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O ' Abstract 

It is pointed out that the dynamics of the order parameter at a 
thermal critical point obeys the precepts of the nonextensive Tsallis 
^ I statistics. We arrive at this conclusion by putting together two well- 

f-*) ■ defined statistical- mechanical developments. The first is that critical 

Tlj" , fluctuations are correctly described by the dynamics of an intermittent 

I nonlinear map. The second is that intermittency in the neighborhood 

' of a tangent bifurcation in such map rigorously obeys nonextensive 

■ statistics. We comment on the implications of this result. 
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The modem theory of critical phenomena has a distinguished history 

|lj of achievements in understanding the source of the scaling and universal- 
ity associated to continuous phase transitions. The introduction of renormal- 
ization group (RG) concepts provided the means for calculating exponents, 
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scaling functions and marginal dimensionalities of critical points in thermal 
and other statistical-mechanical models. Moreover, the RG methods have 
been fruitful in other areas of physics, in condensed matter problems, in 
non- linear dynamics and in other fields |4|. The success of the RG strategy 
in handling problems involving many length scales is illustrated by the use 
of a coarse-grained free energy or effective action. For instance, the equilib- 
rium configurations of Ising spins at the critical temperature shows magnetic 
domains on all size scales and these are suitably studied by means of the 
Landau- Ginzburg- Wilson (LGW) continuous spin model 

Here we may add yet other universal aspect to critical point properties, 
the nonextensivity of order-parameter fluctuations (as explained in more de- 
tail below and in Ref. |5|). This previously unidentified property develops 
as the size of subsystems or domains of a thermal system at its critical point 
is allowed to become infinitely large. These domains have been shown [H]- 
0, via the use of the LGW free energy, to possess intermittency properties, 
and these, in turn, are seen to comply ^21 with the presumptions of the 
nonextensive generalization ^U] QI] of the Boltzmann-Gibbs (BG) statistical 
mechanics. 

In a series of papers [01-0 a bridge has been built connecting the equilib- 
rium dynamics of fluctuations at an ordinary thermal critical point with the 
intermittent dynamics of critical nonlinear maps at a tangent bifurcation. 
With the initial aim of investigating the origin of the relationship between 
the fractal geometry of clusters of the order parameter and the critical expo- 
nents of the phase transition, a theoretical approach was devised that lead 
to the evaluation of the partition function of a cluster or domain where the 
dominant contributions arise from a singularity (similar to an instanton) lo- 
cated in the space outside it |7j . Subsequently [El , a nonlinear map for 
the average order parameter was constructed whose dynamics reproduce the 
averages of the thermal critical properties. This map has as a main feature 
the tangent bifurcation and so time evolution is intermittent. 

Recently we have shown that the known exact static solution 

of the RG equations for the tangent bifurcation in nonlinear maps also de- 
scribes the dynamics of iterates. The fixed-point expressions have the specific 
form that corresponds to the temporal evolution of ensembles of iterates pre- 
scribed by the nonextensive formalism. The proof rests on the derivation of 
the sensitivity to initial conditions exclusively from RG procedures with- 
out approximations followed by comparison with the nonextensive ^j. The 
study of the intermittency transition has been expanded via detailed deriva- 
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tion of their g-generalized Lyapunov coefficients Xq and interpretation of the 
different types of sensitivity [13^. Likewise, the properties of the intricate 
trajectories at the edge of chaos in unimodal maps have been analytically 
obtained leading too to the determination of Xq and interpretation of the 
dynamics at the strange attractor jl5,. 

We combine the results mentioned in the previous two paragraphs to 
point out the nonextensive nature of critical fluctuations. This conclusion 
relates to the quest of the physical circumstances for which BG statistics 
fails to be applicable and its nonextensive generalization might offer correct 
descriptions. These anomalous situations are signalled by the vanishing of the 
Lyapunov coefficients (a single coefficient Ai = for one-dimensional maps) 
and exhibit nonergodicity or unusual phase space mixing P2]-C11 CH]- At 
an intermittency transition hindered or incomplete mixing in phase space 
arises from the special 'tangency' shape of the map at its origin. This has 
the effect of confining or expelling trajectories causing irregular phase-space 
sampling, in contrast to the thorough coverage in generic states with Ai > 0. 
The occurrence of anomalous nonextensive states appears to be related with 
a nonuniform convergence of limits, such as the thermodynamic and infinite 
time limits. Here we comment on the nonuniform convergence associated to 
the description of critical domains as obtained from the LGW free energy. 

Below we expand our arguments. 



2 Criticality and intermittency 

We briefly recall basic elements and results of the approach in Refs. ini-IH]- 
The starting point is the partition function of the (i- dimensional system at 
criticality, 

Z= f D[<p]expi-r,[<p]), (1) 



where 

(2) 



TcM =91 f dV 
Jn 



is the critical LGW free energy of a system of d- dimensional volume fl, cj) is 
the order parameter (e.g. magnetization per unit volume) and 6 is the critical 
isotherm exponent. The partition function Z was evaluated for a subsystem 
of size V and by approximating the path integral in Eq. (H)) as a 

summation over the saddle-point configurations of Fcf^] - an approximation 
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valid for gi ^ 1. Integration of the Euler-Lagrange equation associated 
to the saddle points of Tcf^], and identification of dominant saddle points, 
leads to power-law magnetization profiles for (spherically symmetric) critical 
clusters 0(r), and, finally, to the evaluation of the free energy Fcf^] and the 
partition function Z in closed form [H] [ZI- In doing this it was important to 
notice that only configurations with tq ^ R, where vq is a system-dependent 
reference position tq = ro{g2,S, (f)(0)), have a nonvanishing contribution to 
the path integration jB] |Z|- These configurations vanish for the infinite size 
system, so there is nonuniform convergence in relation to the limits i? — > oo 
and ro — >■ oo, a feature that is significant for our connection with g-statistics. 
Based on the above results the fractal geometry of the critical clusters was 
determined and its relationship with the exponent S was derived. The power- 
law dependence of the magnetization on the cluster radius was identified 
with the fractal dimension of the geometry of the cluster [H] [7j. Multifractal 
properties appear when global, rather than a single cluster, properties are 
considered. 

Subsequent to this development, a link was revealed jS] between the 
fluctuation properties of a critical system described by Eq. ^ and the 
dynamics of marginally chaotic intermittent maps. By considering the space- 
averaged magnetization $ = /y (j){x)dV, the statistical weight 

p($)=exp(-r,[<|.])/Z, (3) 

where rc[$] ~ gig2^^^^ and Z = J exp{—T d^]) , was seen to be the invari- 
ant density of a statistically equivalent one-dimensional map. The functional 
form of this map was obtained as the solution of an inverse Frobenius-Perron 
problem jS] ^Tj. For small values of $ the map has the form 

<l>„+i = $„, + M<l>^+i + e, (4) 

where the amplitude u depends on gi, g2a.iad 6, and the shift parameter 
e ~ R~'^. Eq. ^ can be recognized as that describing the intermittency 
route to chaos in the vicinity of a tangent bifurcation ^H]- The complete 
form of the map displays a superexponentially decreasing region that takes 
back the iterate close to the origin in one step. Thus the parameters of 
the thermal system determine the dynamics of the map. Averages made 
of order-parameter critical configurations are equivalent to iteration time 
averages along the trajectories of the map close to the tangent bifurcation. 
The mean number of iterations in the laminar region was seen to be related 
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to the mean magnetization within a critical cluster of radius R. There is a 
corresponding power law dependence of the duration of the laminar region 
on the shift parameter e of the map [H]. For e > the (small) Lyapunov 
exponent is simply related to the critical exponent 6 0. 

3 Intermittency and nonextensivity 

Next, in a few words, we recall the nonextensive nature of the nonlinear 
dynamics at the tangent bifurcation J2] jHI- At this transition, the inter- 
mittency route to chaos, the ordinary Lyapunov exponent Ai vanishes and 
the sensitivity to initial conditions = \dxt/dxo\ (where Xt is the orbit po- 
sition at time t given the initial position xq at time t = 0) is no longer 
given by the BG law = exp(Ait) but acquires either a power or a super- 
exponential law [13j. The nonextensive formalism indicates that C,t is given 
by the g-exponential expression, 

= exp^(A,t) = [1 - (g - l)X,t]''/^'^''\ (5) 

containing the entropic index q and the g-generalized Lyapunov coefficient 
Xq. Also, according to the generahzed theory, the g-Pesin identity Kg = 
Xq replaces the ordinary Pesin identity Ki = Ai, Ai > 0, where Kq = 
[Sq{0) — Sq(t)] and Ki = t^^[Si{t) — 5*1(0)] are the rates of entropy incre- 
ment based on the Tsallis entropy 

(where In^?/ = (y^^'^ — 1)/(1 — g) is the inverse ofexp^^y)) and the BG entropy 

w 

1=1 

respectively. (Above, Pi{t) is the distribution obtained from the relative fre- 
quencies with which the positions of an ensemble of trajectories occur within 
cells i = 1,...,W at iteration time t). In the limit g — ^ 1 the expressions 
for the nonextensive theory reduce to the ordinary BG expressions. See Ref. 
[T^ for a more rigorous description of the Pesin identity and related issues. 

Assisted by the known RG treatment for the tangent bifurcation JH] , the 
formula for has been rigorously derived ^2] ^H] and found to comply with 
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Eq. (0). Also the validity of the g-Pesin identity has been substantiated for 
this problem [T2] • The tangent bifurcation is usually studied by means of the 
map 

f{x) = e + x + u \x\' + 0{\x\"), u>0, (8) 

in the limit e — 0. The associated RG fixed-point map x' = f*{x) was found 
to be 

x' = x exp,{ux'-^) = x[l -{z- l)ux'-Y^/^'-^\ e = 0, (9) 

as it satisfies = a~^f*{ax) with a = 2^/*^^^^-' and has a power-series 

expansion in x that coincides with Eq. (jSJ in the two lowest-order terms. 
(Above x^~^ = sgn(x)). The long time dynamics is readily derived 

from the static solution Eq. Q, one obtains 

Uxo) = [I - {z - l)axrH]-^/^^-'\ (10) 

and so, g = 2 - and Ag(so) = zaxQ~^ jI2I [121 • When g > 1 the left- 
hand side {x < 0) of the tangent bifurcation map Eq. (jH)) exhibits a weak 
insensitivity to initial conditions, i.e. power-law convergence of orbits. How- 
ever at the right-hand side {x > 0) of the bifurcation the argument of the 
g-exponential becomes positive and this results in a 'super-strong' sensitivity 
to initial conditions, i.e. a sensitivity that is faster than exponential |13j . 

4 Nonextensivity and criticality 

The implications of joining the results described in the previous two sections 
are apparent. In the subsystem of infinite size R ^ oo the dynamics of 
critical fluctuations obey the nonextensive statistics. This is expressed via 
the time series of the average order parameter $„, i.e. trajectories with 
close initial values separate in a superexponential fashion according to Eq. 
ffTUj) with g = (2(5 + l)/(5 + 1) > 1 and with a g-Lyapunov coefficient 
determined by the system parameter values 6, gi, g2 and $o Also, when 
considering an ensemble of trajectories {$n} with prescribed distribution of 
initial conditions, the g-Pesin identity Kg = \q holds with the rate of entropy 
production Kg evaluated according to the nonextensive entropy Sg. 

It is interesting to comment on the conditions for the incidence of g- 
statistical properties at criticality and the manner in which these develop. 
The order parameter profile for a large but finite-size domain i? >> 1 has 
the form [n]-[H] 

0(r) = A,(r2-r2)(2-'^)/2,d>3, (11) 
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that, parenthetically, can be rewritten in terms of a g-exponential. There is 
a singularity in 0(r) when r = tq, the reference position. We keep in mind 
the requirement i? <C tq for the critical clusters 0(r) to be of relevance to 
the partition function Z and also the map shift parameter dependence on 
the domain size e ~ R~'^. The time evolution of $ displays laminar episodes 
of duration < n >~ e~<^/('^+i) and the Lyapunov coefficient in this regime 
is Ai ~ e |7]. Within the first laminar episode the dynamical evolution 
of $ obeys g-statistics, but for very large times the occurrence of many 
different laminar episodes leads to an increasingly chaotic orbit consistent 
with the small Ai > and BG statistics is recovered. As R increases {R -C ro 
always) the time duration of the nonextensive regime increases and in the 
limit R —* oo there is only one infinitely long laminar nonextensive episode 
with Ai = and with no crossover to BG statistics. On the other hand when 
R > tq the clusters 0(r) are no longer dominant, for the infinite subsystem 
R oo their contribution to Z vanishes [n]-[H] and no departure from BG 
statistics is expected to occur. 

This study is developed in Ref. jS]. 

5 Current questions on (/-statistics 

What is the significance of the connections we have presented? Are there 
other connections between critical phenomena and transitions to chaos? Are 
all critical states - infinite correlation length with vanishing Lyapunov co- 
efficients - outside BG statistics? When does BG statistics stop working? 
Where, and in that case why, does Tsallis statistics apply? Is ergodicity fail- 
ure the basic playground for applicability of generalized statistics? Is there 
a link between critical dynamics and glassy dynamics? 

Suggestive results and partial answers to these questions are given here 
and in the cited references. In relation to this it is pertinent to mention 
the following recent development. In Ref. (TO] it is argued that the dy- 
namics at the noise-perturbed edge of chaos in logistic maps is analogous to 
that observed in supercooled liquids close to vitrification. That is, the three 
major features of glassy dynamics in structural glass formers, two-step relax- 
ation, aging, and a relationship between relaxation time and configurational 
entropy, are displayed by orbits with Ai = 0. Time evolution is obtained 
from the Feigenbaum RG transformation, it is expressed analytically via q- 
exponentials, and described in terms of nonextensive statistics. 
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